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ABSTRACT- * Sardina pikhardus populations from the Aegean and Ionian seas have been 
examined for differences in growth rate using the statistical test described by Kappenman 
(19SI). Results suggest that growth in both populations can be described by the same von 
Bertatanffy growth curve. Further studies of the spawning populations are needed to explain 
the one year difference observed in recruitment time between these populations, and to clarify 
whether they constitute different stocks. 


RESUME- - Le lest siatistique deed! par Kappenman (19$ I) a etc utilise pour etudier les 
differences entre les taux de crotssance des populations de Sardinia pikhardus des mers Egee 
et lonienne, Les resuliats mettent en evidence que la croissance de ces deux populations peut 
etre decritc par la me me courbe de von BertalanfFy, Des etudes ulteneures sur la reproduc¬ 
tion seraient necessaires pour expliquer la difference d un an au niveau du recrutement des 
deux populations et pour savoir s i] s agit de deux stocks differents. 
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The sardine (Sardina pilchanius Walb„ 1792) is one of the most exploited 
pelagic fish in the Greek seas. Although it inhabits the coastal and offshore waters 
in the entire area of Lhc Aegean and Ionian seas, the principal fishing grounds are 
situated in die northern Aegean, while the main fishing grounds in the Ionian sea 
are located in die central part of the area. I he fishes in both areas are caught by 
identical purse-seine techniques. 

Monthly random samples of sardine were taken during the fishing seasons 
(March-October) of 1983 and 1984 from commercial catches unloaded at the main 
quays of the fishing ports of Kavala and Patra located in the northern Aegean and 
central Ionian seas respectively. They showed a noticeable difference in their length 
composition, such dial catches from Kavala consisted of larger individuals than 
those from Patra in almost all samples taken (Figs L 2), although the fishing 
method is identical in both areas. 

It is known that within a unit stock there can be large and consistent 
differences in size and age com position between catches from different fishing 
grounds. The genetic findings by Spanakis et al. (1989) that sardines from the 
Aegean and Ionian seas do not form one panmictic population, in conjunction 
with the above mentioned observations of their length composition, stimulated us 
to examine whether there are differences in growth rates between the sardine popu¬ 
lations of the two areas. 
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Fig. 1, - Length Frequency distribution of sardine samples taken from the ports of Patra 
(continuous line) and Kavala (dotted line) in 1983. 

Several attempts have been made in the past to provide information on 
sardine growth in the Mediterranean (Murat, 1935; Rodriguez-Roda and 
Larraneta, 1954, 1955; Sinovcic, 1984: Bouchereau and Ghazi, 1988); however, no 
growth studies exist for the populations in the Aegean and Ionian seas. This paper 
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Fig. 2. - Length frequency distribution of sardine samples taken from the ports of Patra 
(continuous line) and Kavala (dotted line) in 1984, 


is primarily concerned with the evaluation of the observed differences in length 
composition of sardine populations from the above areas, in terms of growth rate, 
but also w ith the actual growth of sardines in these areas. 
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MATERIAL AND METHODS 

A thousand sardines were obtained from commercial purse-seine catches 
from the fishing ports of Patra and Kavala in the Ionian and Aegean seas. 
Samples were taken monthly (50 individuals from each port) from March 1983 to 
October 1984. Animals were taken randomly from the catches unloaded in the 
quays previously to their separation by stee. Lor each fish the total length (TL) to 
the nearest mm was measured, the sex was determined by gross examination of the 
gonads, and several scales were removed from the region below the pectoral fm. 
Since sardines loose their scales very easily, if scales were not found in dial region 
then others were taken from nearby. Six scales from each specimen were cleaned 
with water, placed on a 0 J mm thick cellulose acetate plastic plate and pressed at 
a temperature of 95 S C for two minutes. In this way the prints of die scales were left 
on the plate when die scales were removed. The prepared plates with die prints of 
the scales were read using an image processing system, consisting of a computer 
equipped with a frame grabber, a video-camera, a monitor and a software package 
(Anonymous, 1987), The image of the scales was captured by die video camera, 
digitized by the frame grabber and displayed on the monitor. It was possible to age 
the scales automatically because the software was able to judge the different 
luminance values of the various image areas and then mark certain bright or dark 
zones. 1 lowever, in practice, it was necessary to interfere manually in order to 
eliminate false marks or to add missing ones. All prints w ere read in order to find 
the most readable one which was used for age determination. The number of 
annuli from this print was used to estimate die age of each specimen. Scale prints 
were read twice by a single reader. To reduce subjective bias in age estimations, 
where possible, an interval of 3 months separated the two readings. When the two 
readings of the same print resulted in different age estimates, the scales were 
considered to be unreadable. 

Sex composition of the samples was compared for independence of 
sampling port using the chi-square test. 

The von Bertalanffy growth curves for the two supposed populations w'ere 
compared using the method described by Kappenman (1981), i,e, the data sets 
from the two areas were used in order to decide which of two models, M ( or \1 2 , is 
preferred; in Mj the growth curves for the two populations are the same, while in 
M 2 the two populations have different growth curves. According to this method, for 
each particular age-length observation Lwo curves were constructed (one for each 
model), without taking into account this observation, and the predicted lengths in 
each model, as well as the squared differences between predicted and observed 
lengths have been computed. The latter have been summed and the model giving 
the smallest sum was adopted. In the present study the mean lengths aL age w'ere 
used instead of the raw age-length data. However this modification does not affect 
the reliability of the Kappenman method. Differences in mean observed lengths at 
age for the samples from the different ports were tested for significance using single 
classification analysis of variance (ANOVA), Finally, estimates of theoretical 
growth in length under the preferred model were obtained by fitting the observed 
mean length at age data to the von Bertalanffy growth equation using the Newton 
iteration method (Sparre, 1987), All statistical inferences were based on a 
significance level of a = 0.05, 


RESULTS 

The sex ratios for all sampled individuals and for the ones which had been 
aged are shown in table L In both cases the chi-square test detected a slightly 
significant difference in sex composition of samples taken from the two ports (x a = 
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4.97 and 4.09 respectively, where x 2 M — 3.84). Therefore, specimens of both sexes 
were combined for age and growth estimations* 

Table 1. - Sex ratios of sampled and aged fish. 


AHEA 

FEMALES 

MALES 

XAVALA 

PATMA 

sampled 

Aged 

Sampled 

Aged 

297 (59.4*) 
262 {52.41} 

163 {54*3%) 
119 {45.31) 

2S3 {40.6%) 

233 (47,61) 

137 {45,7%) 
141 (54*21) 


The calculated length at age obtained under models Mj and M 2p compared 
with the corresponding mean observed lengths (Tables I1 T 111) by using the 
Kappcnman test, showed that growth of the populations from both seas can be 
described hy the same curve since the sum of squares of the difference between 
observed and predicted lengths under model M, was found to be smaller than the 
one under \1 2 (2Q21.94 and 2820.96 respectively)- The growth curve for both po¬ 
pulations is given by the following equation: 

L x = ISI*23[l-e^ l + ^ 

where L t = length at age t* 

Table IV shows the mean observed lengths at age for both areas and the 
ones predicted by the common von Bcrtalanffy growth equation. Table IV also 
shows that die samples from Kavala pooled for the years 1983 arid 1984 consisted 
of older individuals than those from Patra. It was found that the mean observed 
age of sardine samples from die Aegean sea was 2.85 years while that from the 
Ionian was 1.87* 

The above mentioned difference in age composition occurs not only for the 
pooled data hut also for those in each of the sampling years 1983 and 1984. 


Table II. - Ages, mean observed total lengths and predicted lengths (mm), under models Ml 
and M2, for samples from Patra. 


AGE 

MZATt CSS* 71 

9RECIC71S 71 
UltSER M 2 

?3£5XCT£X> TL 

3*23 

120-53 

117 . 3 a 

114,54 

0*42 

122.00 

117.55 

iia.91 

♦ Q. 67 

119 * 50 

134* 55 

125.32 

0,75 

115.CO 

129 *5S 

127*76 

1*53 

134*00 

152.57 

154.75 

1*42 

139*73 

130.74 

156*52 

1.50 

145 * 50 

131.17 

136*31 

1, EE 

141.19 

135*36 

138,95 

1*57 

150,62 

158.54 

140*19 

1.75 

123,00 

133,43 

141-29 

1*33 

146,50 

137,54 

141-21 

2*33 

143 *57 

144.76 

147.08 

2*42 

147.72 

145,24 

147*57 

2*50 

153 * 14 

143.53 

143,59 

2-=a 

146.50 

143*22 

149,32 

2*75 

144.CO 

150.29 

>zS 

5,33 

149,57 

146.73 

152.49 

5- 55 

147.SO 

154.25 

135♦46 

3,42 

153*33 

155*49 

155*50 

3* =3 

153,36 

1=3*14 

157,06 

5*55 

155,30 

153,14 

137,61 

3 * 7 ; 

143-30 

151*58 

*33,40 

5 * 32 

255,00 

161*37 

150,13 

4,42 

155*30 

154* £5 

- j * _3 

4*52 

172.SC 

151,S5 
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Table 111, - Ages, mean observed total lengths and predicted lengths (mm) under models Mj 
and M 2 , for samples from Kavala. 


AGE 

MEAN OSS. 7£ 

PREDICTED T! 

UNDER 

PREDICTED ”L 

UNDER 

I. IS 

140.30 

140.07 

133*51 

1*33 

143*40 

138.18 

136*24 


138-46 

134.51 

136.39 

1*50 

140*71 

144.32 

137.94 

1*67 

144*00 

143*37 

133* 71 

1*75 

142*00 

143.17 

140*30 

1.33 

139*00 

145.11 

142.23 

2,08 

156.15 

145*10 

143,68 

2*25 

145.63 

148*99 

145.74 

2*33 

149-79 

143.77 

146.10 

2*42 

150*82 

152.02 

147,15 

2,SO 

151-42 

153.30 

147.79 

2,58 

153.33 

145*63 

149,08 

2,67 

147.65 

153*27 

150,86 

2 * 75 

149.35 

153*13 

161.97 

2,83 

144,17 

153-13 

152,43 

j.oe 

162 * 13 

162 * 65 

154*10 

3.25 

154*33 

158,21 

154.37 

3.33 

157,17 

153.46 

155.96 

3.42 

159.50 

169.12 

156-60 

3*50 

164*52 

155*65 

156*99 

3.55 

153*30 

160*99 

157.53 

3*57 

156.33 

161,07 

153.40 

3*75 

155.00 

153-59 

159*23 

3.33 

153.60 

158.54 

160*39 

4* OS 

170.00 

153.46 

160-76 

4.33 

170,00 

164*29 

162.64 

4.42 

166-00 

160*35 

162,58 

4-50 

162*00 

164,43 

162.E6 

4,53 

165.00 

162 - 50 

164*34 

4.75 

164.25 

163-49 

165*09 

4*33 

161,00 

103,91 

156 -4C 


Table IV. - Age group frequencies and mean total lengths observed and predicted by the von 
Bertalanffy length at age for aged fish. 


AREA 

ACS 

FREQUENCY 

085* TE 

BERT. TE 



1983 

1984 

Total 




0 

36 


40 (15,411 

119-27 

112*05 


1 

40 

72 

112 (43.1%| 

137*23 

129,98 

Panra 

2 

33 

44 

82 (31*5%) 

147.38 

143.26 


3 

16 

7 

23 ( 8.3%) 

154.61 

153.10 


4 

3 

0 

3 { 1,2%) 

163*50 

160*39 


0 

0 

□ 

0 

- 

112.05 


1 

32 

11 

43 (14*5%) 

141.20 

129*93 

Ravala 

2 

60 

74 

124 (44.7%) 

149.31 

143.26 


3 

49 

43 

98 (32*7%) 

156*76 

153.10 


4 

6 

19 

25 ( 8,3%) 

165.46 

160.39 


Specifically* it was found Lhat the mean observed ages of samples obtained in 19S3 
from the ports of Kavala and Paira were 2*70 and LS2 years respectively, while 
the corresponding ages of samples obtained in 1984 were 3.00 and 1.93 years. 

Considering mean observed length at age between samples from the tw^o 
areas A NOVA delected no significant difference after age 1 (Table V)* 
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AG3 

F-RATIO 

CRITICAL F-VALUE 

X 

20,53 

3.91 

2 

0* 57 

3,39 

3 

0- 63 , 

3.92 

4 

0*06 

4,22 


Tabic V, - Results of ANOVA tests performed 
to detect differences in mean observed lengths 
at age between samples from the different 
ports. 


DISCUSSION 

Several authors have pointed out that estimation of the von BertalanfTy 
growth parameters is highly affected by the composition of the sample (Paloheimo 
and Dickie, 1965; Southward and Chapman, 1965; Pauly, 1979). Therefore, 
comparison of curves obtained from samples with different length compositions 
makes necessary the application of a proper statistical method. The Kappcnman 
lest used in the present study is simple in its concept and seems to give reasonable 
results. It has shown that growth in both populations studied is represented by the 
same von BertalanfTy growth curve. The back calculation method was not applied 
in the present study since scales were not taken from the same body area in all 
specimens, thus they differ in size and shape. 

The use of an image analysts system for scale reading was found to be more 
convenient, faster and gave more accurate measurements of scale and annulus 
radius in comparison w r ith the traditional w T ay of reading under a microscope. 
Moreover, the ability of the present system to store the results made the subsequent 
work of data processing easier. 

Considering the differences found in length and age composition of catches 
from the two areas, table IV shows that they may be due to the fact that 
recruitment of sardines in the central Ionian sea starts one year earlier than in the 
northern Aegean sea. However, since the fishing season is the same in both areas 
and fishing is carried out by using the same gear, the exploitation pattern should 
not affect the time of recruitment. 

Certainly a difference in time of recruitment is not unusual for populations 
distributed over a wide geographical range, but the difference of almost one year 
observed in the present study may need some additional explanation. 

Spanakis et aL (1989) state that though their findings suggest that sardine 
populations from the Aegean and Ionian seas do not form a panmictic population, 
this is not to say that these populations constitute discrete stocks. Thus their 
findings and the results of the present study arc in agreement with Gulland's (1983) 
view that average age or size is a poor parameter to use for stock differentiation. 

Obviously further studies on the spawning populations are needed to clarify 
whether these populations do or do not constitute different stocks. 


Table VI. - Length (mm) at age from different studies in the Mediterranean. 


AGE 

Fage 

(1920) 

(Goife 

du 

Lion, 

France} 

Murat 
(1935) 
(Caati* 
-gliaHe, 
Algeria} 

Rodriguez 

and 

Larrarteta 

(1954) 

(Castfiilon, 

Spain) 

Sinovcic 

(1934) 

(Central 

Adriatic) 

Boucherau 
and Chazi 
(1909) 
(Algeria) 

Present 1 

study 

t 

98 

110 

141 

102,13 

102.93 

129.93 

2 

129 

145 

155 

140.09 

123,39 

143,25 

3 

141 

1S5 

139 

154.02 

140.S3 | 

153.10 

- 

149 

150 

194 

179-13 

153,99 

ISO.39 
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Concerning growth rale Usclf, table VI summarizes results from different 
growth studies on sardine in the Mediterranean. Estimated values of length at age 
from the present study lie between the corresponding values reported for other 
Mediterranean areas. However, die scarcity of small individuals in the present 
study may have affected the estimation of the von Rcrtalanfly growth parameters 
and resulted in overestimating length, mainly at age 1, but this does not affect the 
method applied for growth curve comparisons (Kappenman, 1981). 
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